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Abstract 

In this letter we present the three-loop virtual corrections to the Higgs boson 
production in the gluon fusion channel where finite top quark mass effects are taken 
into account. We perform an asymptotic expansion and manage to evaluate five 
terms in the expansion parameter Mfj /Mf . A good convergence is observed almost 
until M H « 2M t . 
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1 Introduction 



Among the main tasks of the CERN Large Hadron Collider (LHC) will be the uncovering 
of the mechanism which provides particles with their masses. A crucial role in this respect 
is assigned to the Higgs boson whose discovery is awaited with great eagerness. 

At LHC the Standard Model Higgs boson is mainly produced in the so-called gluon fusion 
process where two gluons couple via a closed quark loop to the Higgs boson. The leading 
order (LO) process for this channel has been evaluated in Refs. [1-4] and already almost 
15 years ago also the next-to-leading order (NLO) corrections became available [5,6]. To 
this order the production cross section could be evaluated without any assumption on 
the hierarchy between the mass of the quark in the loop, the Higgs boson mass, and the 
partonic center-of-mass energy. 

On the contrary, at next-to-next-to-leading order (NNLO) only quantum corrections in- 
volving the top quark Yukawa coupling are available. It has been performed under the 
assumption that the top quark is much heavier than the Higgs boson. In this limit it is 
suggestive to construct an effective theory where the top quark is integrated out. The 
coefficient function of the corresponding effective operator has been computed to NNLO 
in Refs. [7,8] (see also Ref. [9]) and the production cross section has been evaluated in 
Refs. [10-13]. Let us mention that recently the virtual contributions to the NNNLO 
corrections have been completed [7, 14] . 

In this letter we provide the first results beyond the effective-theory approach where 
three building blocks are required to NNLO: virtual three- loop corrections to the 2 — ► 1 
process, two-loop corrections to the 2^2 process where next to the Higgs boson a 
parton is radiated off, and one-loop corrections with radiation of two additional partons. 
We present results for the three-loop virtual corrections to the process gg — > H including 
finite top quark mass effects. Our results constitute a building block for the NNLO 
corrections beyond the heavy top quark limit. 

Let us mention that there are also results beyond the fixed-order approximation. In par- 
ticular, in Ref. [15] large logarithms in connection with soft gluon radiation have been 
resummed. A step further has been taken in Ref. [16], where certain 7r 2 terms have been 
resummed leading to a perturbative series which is significantly better behaved as com- 
pared to the unresummed approach. In Ref. [17] the limit of high partonic center-of-mass 
energies has been considered for the gluon-fusion process and an approximation for the 
NNLO cross section has been derived which goes beyond the M t — > oo result. Electroweak 
corrections have been considered in Refs. [18,19]. For recent numerical predictions of Higgs 
boson production in gluon fusion both at the Tevatron and the LHC we refer to Ref. [20]. 

The remainder of the paper is organized as follows: in the next section we briefly describe 
details of our calculation and the Section [3] contains our results and conclusions. 
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Figure 1: Sample diagrams contributing to the NNLO virtual corrections to gg — > h. 

2 Calculation 



The existing NNLO calculations to the Higgs boson production have been performed in 
the framework of an effective theory where the top quark has been integrated out. In 
this way effective vertices^ are generated between the Higgs boson and two, three or four 
gluons. The number of loops to be considered for the calculation of the cross section 
is effectively reduced by one, leading to virtual two-loop 2 — > 1, one- loop 2^2 and 
tree-level 2 —>■ 3 corrections. 

In contrast to this approach, our starting point is the full QCD with six active flavours. 
Some sample diagrams contributing to the virtual corrections are shown in Fig. [TJ al- 
together 657 three-loop diagrams have to be considered which in the sum lead to the 
structure 



where p = Mjj/M^ and q\ and q2 are the incoming momenta of the two gluons with 
polarization vectors £^{qi) and e u (q2)- In our calculation we construct projectors on X 
and Y and check the condition X = — Y which follows from gauge invariance. The 
ellipses in Eq. (JT|) represent further structures which do not contribute to the physical 
cross section. They would receive contributions from vertex corrections with external 
ghosts which we do not consider in this paper. 



The virtual contribution to the partonic cross section can be cast in the form 



a 



virt 



0"LO 



(l+^^) + (^) 2 ^)+...) , 



(2) 



where the LO cross section is given by 

Gf a 2 s /o0,e) 



o"lo 



6{1 -x), 



(3) 



288\/2V (1 - 

with x = Mft/s and dimension of space D = 4 — 2e. is the partonic center-of-mass 



1 Thc effective coupling has even been computed to four- and five-loop order in Refs. [7] and [21,22], 
respectively. 
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Figure 2: Factorized regions appearing from the asymptotic expansion of the double-scale 
integral (a). Solid lines carry the mass M t , dotted lines mass Mh, and the dashed lines 
are massless. Cases (b),(c), and (d) correspond to one, two, and three loop momenta at 
the scale M t , respectively, with the remaining loops at the scale Mjj. 
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where an expansion through 0(p 4 ) and 0(e 4 ) has been performed in the second line. In 
Section [3] we will present results for d^ 1 ' and . 

Note that the quantities 5^ only depend on p, since the Feynman diagrams have to 
be evaluated for on-shell external partons. It is thus tempting to evaluate them in the 
limit Mh <C 2M t which is expected to show good convergence properties even up to 
M H m 2M t [23]. In Ref. [24] the NNLO corrections to the decay of a Higgs boson into 
gluons have been considered. The optical theorem in combination with the asymptotic 
expansion was used in order to evaluate three expansion terms in M^/M 2 where rapid 
convergence has been observed for M# m M t . The asymptotic expansion [25] in the limit 
Mh 2M t leads to one-, two- and three-loop vacuum integrals where the scale is given 
by the top-quark mass and to one- and two-loop vertex diagrams with massless internal 
lines and external momentum at the scale Mh- In Fig. [2] we exemplify the asymptotic 
expansion in diagrammatic form for a typical contribution. 

For our calculation we have used two independent set-ups. In the first one all Feynman 
diagrams are generated with QGRAF [26]. The various diagram topologies are identified 
and transformed to FORM [27] notation with the help of q2e and exp [28,29]. The program 
exp is also used in order to apply the asymptotic expansion (see, e.g., Ref. [25]) in the 
various mass hierarchies. The actual evaluation of the integrals is performed with the 
packages MATAD [30], which is used for the vacuum integrals, and FIRE [31], employed to 
reduce the massless three-point functions to master integrals. The latter can, e.g., be 
found in Ref. [32]. 

The second set-up also relies on QGRAF for the generation of the Feynman diagrams. 
Afterwards the asymptotic expansion is done with a Perl program, and, two- and three- 
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loop integrals are reduced by an independent implementation of the Laporta algorithm [33, 
34]. 

We have performed the evaluation of the vertex corrections up to order p 2 for general 
QCD gauge parameter and checked that it drops out in the sum of all bare three-loop 
diagrams which serves as a welcome check of our calculation. 

In the sum of all three-loop diagrams we observe poles up to order 1/e 4 which is due to 
a mixture of ultra-violet and infra-red singularities. As usual, the ultra-violet poles are 
treated via renormalization. In our case we have to renormalize the gluon wave function, 
top quark mass and strong coupling constant to two-loop order. The remaining infra-red 
poles are only cancelled after including the corrections from the real radiation and mass 
factorization. 

(5) 

In the following section we present our results expressed in terms of a s , the strong 
coupling constant in the MS scheme defined in five-flavour QCD, and the on-shell top 
quark mass. Since our two-loop result contains poles up to 0(1/ e 2 ), the one-loop top 
quark mass counterterm is needed up to 0(e 2 ) which can be found in Ref. [35]. The 
two- loop counterterm for a s can be found, e.g., in Ref. [36] and the one for the gluon 
wave function in Ref. [7]. 

The transition from to ai 5 "* is performed with the help of the formulae derived in 
Ref. [7]. Since there are poles in the NLO expression, higher order terms in e are necessary 
for the decoupling relation. The explicit result can be found in Eq. (12) of Ref. [37]. 



3 Results 

At the three-loop order we were able to evaluate the first five terms in the expansion 
around p = 0. The NLO and NNLO corrections to the partonic cross section are given 
by (adopting common MS conventions) 
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Figure 3: Finite part of 5^ (left) and 5^ (right) as a function of p. The longer-dashed 
lines include successively higher orders in p. 
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where L^h = ln(p 2 /M H ), L^t = ^ n {M H /M 2 ) and ((n) is the Riemann's zeta function. 
Furthermore, SU (3) colour factors have been applied and the number of massless quark 
flavours is set to rii = 5. The analytic expression for the generic values of N c and rii can 
be found in [38]. 

We have checked that the expansion of 5^ to order p agrees with [39]. At the NNLO the 
leading term in the inverse top quark mass expansion agrees with the results of Refs. [7,10]. 

Although the final result is divergent, it is instructive to look at the finite parts of S^ 1 ' and 
5^ 2 \ In Fig. [3] the corresponding two- and three- loop expressions are shown for p = Mh in 
the range between p = and p = 4 corresponding to Mh = 2M t . The longer-dashed lines 
include successively higher orders in p up to order p 4 . One observes good convergence up 
to p 3 which corresponds to M# ~ 1.7 M t . 

To conclude, we have presented the virtual corrections to the partonic cross section gg — > 
H including finite top quark mass effects. Our calculation confirms the results obtained in 
the framework of the effective theory and provides four more expansion terms in M H /M 2 . 
We observe a rapid convergence almost up to Mh ph 2M t . The results presented in this 
letter constitute a building block for the NNLO corrections to the Higgs boson production 
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in the gluon fusion channel beyond the heavy top quark limit. 

When this paper was in preparation, we had a chance to learn about the parallel publica- 
tion [40] and establish the full agreement of the results. 
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